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ABSTRACT 

We use standard polynomial expansion technique to show the existence of a relation 
between polytropic model and the description of gas spheres at finite temperature. 
A numerical analysis is made concerning the obtained perturbative parameters. It 
is shown that there is a correspondence between polytropic and gas sphere thermal 
models for fermions and bosons. For instance, the polytropic index n displays evident 
correlation with temperature and chemical potential. 
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1 INTRODUCTION 

As they barely emit light, compact objects are considered as 
part of the dark matter in cosmology. Dark matter is a parcel 
of mass that can only be detected by its gravitational effects. 
With this definition and depending on where the referential 
is placed many things can be considered as dark matter; 
for instance, even Jupiter can be considered as dark matter 
for a distant observer. This subject has increasingly been 
accepted as decisive for understanding the evolution of the 
universe. There are even speculations (jPortilho 2003 ') that 
in the Earth, the Chandler wobble excitation and damping, 
one of the open problems in geophysics, can be considered 
as a consequence of geological interaction with an oblate 
ellipsoid made of dark matter. The same phenomenon is 
also observed in other planets. 

Differently from dark energy, which has no effective 
mass but presents its effects similarly decisive for the evolu- 
tion of the universe, dark matter has mass and may consist 
of many kinds of particles: dust, neutrinos, neutrons, pro- 
tons (hydrogen), and even bosons, like a lpha particles, Higgs 
bosons or axions iJKolb fc Turner 199d ) . Many of those par- 
ticles come from the early universe, encompassing some of 
the lore of its evolution. Fraction of those particles are roving 
over the free space. But, in the early universe, some of those 
particles departed from roaming, clustered and formed self- 
gravitating compact objects, like fermion stars and boson 
stars. 

Fermion stars is a general name denoting particular ones 
like neutron star and white dwarf stars. Since Oppenheimer 
and Volkoff (|OpDenheimer fc Volkoff 1939i ). these compact 
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objects have received large attention and many of their 
properties are already determined, and astronomers can use 
the theoret ical information to detect them. From Chan- 
drasekhar (jChandrasekhar 19391 ) one can also learn that 
neutron stars and white dwarf stars can be studied as a 
degenerate Fermi gas under itself gravitational field, with 
an equation of state to determine its pressure and density 
of energy. 

In c ontrast to fermion star there is the so-called bo- 
son star (IRuffini fc Bonazzola 19691 : iLiddle fc Madsen 1992 : 
iMielke 199 ll ). built up with self-gravitating bosons at zero 
temperature. In this case, one can also use relativistic ap- 
proach, but the star structure can be studied directly from 
the stress-energy tensor since it is possible to write a La- 
grangian for bosons in place of using an equation of state. 
Despite this compact object has not been detected yet, there 
are many theoretical effort s to understand its prope rties, 
like formation and stability (iGleiser 198i: lMielke 20031 '). ro- 
tation (de_Sousa 1995; de Sousa 2006), and even the inter- 
action (|de Sousa 19981 ) between bosons and fermions in the 
same spherical system. 

The aim of the present article is to show that there 
is a relation between the statistical mechanics of the gas in- 
side the star and the polytropic model for both fermions and 
bosons, starting from the isothermal gas sphere, and depart- 
ing from this stage to slightly reach a non-zero temperature 
in first approximation. 



2 SYSTEMS OF SELF-GRAVITATING 
BOSONS AND FERMIONS 

It is well known that neutron stars are actually to be con- 
sidered as fermion stars. Neutron stars are kind of self- 
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gravitating Fermi gas and the study of this kind of quantum 
fluid is in the concern of statistical mechanics. Thus, the 
neutron star matter, which is usually treated in the classical 
regime, requires further studies in the quantum regime to be 
better understood specially under such high pressures and 
curvatures as those experienced in stars. 

In Physics there are two kinds of quantum statistics: 
fermion and boson gases. Statistica l mechanics is con sid- 
ered by Ingrosso and Rufhni ( Ingrosso fc Ruffini igT^ ) in 
the context of both fermion and boson stars. If the fermion 
or the boson gas becomes gravitationally bound and sta- 
ble, the gas will reach a characteristic energy density p and 
pressure p obeying the statistics related to temperature and 
average energy, such that: 



16^25'!/? 



5/2 



i^3/2(e,/3) + 1^5/2 (e,/3) 



(1) 

(2) 



p = 3^S2/3'/'[i^i/2(e,/?)+/?F3/2(e,/?)] 
where: 
e = ^i/kT 
/3 = kBT/mc 

5*2 = (2s + l)mV767r^?i^ 

and, m is the mass and s is the spin of the considered particle 
(fermion or boson). The F functions are m athematical tools 
from statistical mechanics (jPathria 1972h . given by: 
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(4) 



where the sign +(— ) refers to fermions (bosons). The gas 
presents a chemical potential /i, and ks is the Boltzmann 
constant. 



3 POLYTROPES 

The study of polytropic stars provided great simplifications 
concerning the treatment of compact objects. In thermody- 
namics polytropes are paths similar to adiabatics, isobarics 
and isothermals, and families of Lane-Emdem equations can 
be separated concerning their polytrope indices n. This is a 
very important tool to classify self-gravitating objects and 
defining their i nternal energy, gravitationa l energy and pos- 
sible stability ijZeldovich fc Novikov 197l| ). 



3.1 Relativistic stars 

In such to consider gas spheres, we use the spherical metric 
element: 



ds 



~B{r)dT + A{r)dr + rdO^ + r sin^ Odcp'^ 



(5) 



Chandrasekhar proposed that a compact object could 
be treated as a perfect fiuid described by an energy- 
momentum tensor: 



Tfiu = PSm" + (p + pc^) Uf^u^ 



(6) 



where u" is the four-velocity of the gas, p is the pressure 
and p is the energy density. Taking units where c = 1, the 
four- velocity vectors u obey the relation u^u^g^^ = —1, in 
such that Wr = We = M0 = and ut — — (ff**)~^^^ = — \/B. 
The non-vanishing Ricci components, Rrr,Ree and Rtt, give 
rise, respectively, to: 
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It is also useful the equation for the hydrostatic equi- 
librium: 



B^ 
B 



2p' 



p + p 
and the equation for the total mass inside the star: 

cR 



M = 



Aixr p{r)dr 



(10) 



(11) 



So, it is possible to define A(r) using the limit to the infinity: 
1 



A{r) 



2GM 



(12) 



Equation ^ can be rewritten as: 
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2GM , ^ 2 
H 47rG'pr 



= — 47rG' (p — p)r 
This last equation can be written as: 

AivGr'^p 



r'^p'{r) = GMp' 
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2GM 



(13) 



(14) 



This equation can be used to determine pressure evolu- 
tion in r. It can be used in Astrophysics applications within 
the Newtonian approach after some non-relativistic correc- 
tions. Two isentropic special cases are of interest: stars at 
absolute zero (neutron stars, white dwarf stars and boson 
stars), and stars in convective equilibrium. 



3.2 Semi-relativistic approach 

When consider ing the Newtonian case, one can take 
(| Weinberg 19721 '): 



P< P 

47rr^p < M 



2GM 



< 1 



(15) 
(16) 

(17) 



Using these approximations in ean. (|14p simplify to: 
rV = GMp (18) 
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which by means of egn. dlll) gives: 



d_ 
dr 



r I 



'4:TTGr p 



(19) 



This equation could be a linear differential equation (LDE), 
except for the fact that it is to be solved both on p{r) and on 
p(r). But, if p and p have a relation (known as 'equation of 
state') then ean. (|19p becomes an actual LDE, with solution 
depending only upon p or p. In fact, this equation of state 
is the polytropic equation: 



Kp'' 



(20) 



where K and 7 are constants. Now one can see that with the 
adequate boundary conditions, such as p'(0) — (to keep 
p(0) finite), ean. (|19|l gives p = p{r). 

This relation is famous for long time since, for instance, 
W. Thomson (Lord Kelvin) in 1887 studied the natural stir- 
ring produced in a great free fiuid mass like the Sun while 
it is cooling at its surface. At Joule's suggestion, Kelvin 
also studied this natural stirring of a moist atmosphere con- 
densation of vapor in the upward currents of air, which is 
nowadays of recogniz ed importance in met eorology, known 
as polytropic change (|Chandrasekhar 19391 ). Following this 
historical tradition, in Astrophysics any star for which the 
equation of stat e takes the form (|20|l is called a polytrope 
l| Weinberg 19721 '). Some typical cases are: 7 = 6/5 is in the 
range of super large gaseous stars; white dwarfs (and gen- 
erally fermion stars) present 4/3 ^ 7 ^ 5/3, where 7 = 4/3 
correspond to largest mass white dwarfs and 7 = 5/3 to 
small mass white dwarfs; there are also the incompressible 
stars with very high polytropic indexes (7 — >■ 00). 



3.3 Isothermal gas spheres 

Isothermal gas spheres are important in Astrophysics since 
it is a starting point to understand composite stars, and 
also to the study of stars consisting of envelopes with dif- 
ferent temperatures llChavanis 20021 ). For a standard star 
IChandrasekhar 19391) . we have from the theorems of the 
equilibrium of the star: 



pHj^^^r^' 



(21) 



where ks is the Boltzmann constant, p is the mean molecu- 
lar weight, H is the mass of the proton, and a is the Stefan- 
Boltzmann constant. 

If the star is in gravitational equilibrium, one can write 
its equation of state as: 



p = Kp + D 
with: 



(22) 



K = 



pH 



D 



3 



Notice that both K and D depends on the temperature. 
Comparing eqns.(|20p and H22p we also can learn that the 
isothermal sphere case is close to correspond to a polytropic 
equation with 7 = 1, if we take D — >■ 0. 



4 TAYLOR EXPANDING THE EQUATION OF 
STATE 

In this article we assume the star is under hydrostatic equi- 
librium, and using equations ([l}-((2)| one can show that: 



V 
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I6S1 



p + s^s'i/?^'' V5/2 - le^Si/?^/ Vi 



/2 



(23) 



V 3S2 

which is a polytropic-like equation, unless for the two ad- 
ditional terms at the end. It can be compared to equation 
(|20|l if we consider: 



Kp^ + cis,e,p) 



(24) 



where K is constant for a star in equilibrium, and the power 
factor 7 is sometimes expressed using 7 = (1 -I- 1/n) or n = 
(7 — 1)~^, where n is called the polytropic index. 

New physics can be inferred if equations (|20|l and 
are generalities of: 

1+'' 



Kp^ 



(25) 



Some considerations about the results so on. The func- 
tion C depends on temperature, but the polytropic equation 
does not explicitly takes into account the temperature of the 
star. Moreover, the case 7 = 1 corresponds to an isothermal 
sphere at constant temperature and can be considered as 
a special situation in Astrophysics, since n — > 00. Hence, S 
takes into account deviations from this case (if 5 = we 
recover the isothermal sphere). 

We can Taylor expand (|25|) around the central density 
p = po- 

p ^K{pI+' + i5 + l)p^p-{5 + l)pi+'}+0{{p-pof,S^) (26) 

Thus, up to first order in p — po and 5: 

p = [K{S + l)pf,] p + Kp]+' - K{5 + l)p]+' (27) 

A comparison between this last equation and (|23|) 
brings us to a system of non-linear equations, which is 
straightforward solved giving: 



2a -1 

3x/25'2/3='/'fi/2 
8%/25i/3^/^F5/2 

(3x/252/33/2Fi/2)' 



where: 






(28) 



(29) 



The value of 5 brings relevant information, since one can 
associate it with the polytropic index, i.e., 5 = 1/n. Apart of 
-Ffc functions one can write: a ~ m?c^ /k^T^ . Thence, since 
function 5 is directly related to a, which depends on the 
mass of the particle on consideration, results for bosons and 
fermions may present large numerical differences. 



5 RESULTS 

The aim here is to show that there is a relation between 
the polytropic model and the statistical mechanics of the 
gas inside the star. Consequently, in the previous section it 
has been shown that the parameters 5, po and K, can be 
related to the temperature and to the chemical potential. 
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Figure 1 . Plot of 5 in function of the temperature T and ciiemical 
potential fj. for fermions. Notice tliat the values of 5 are in order 
of lO^'^, and decreases strongly for higher temperatures. 



Figure 2. Plot of S for bosons. Notice that the values are in order 
10'^^, and decreases as temperature increases. 



We separate the results for bosons and for fermions in such 
to compare the behaviour for these two cases. 

To illustrate the parameters behaviour we have used 
neutrons as fermions, and the Z° bosons. Their masses, spins 
and ot her constants can b e recovered, e.g., in Particle Data 
Group (|Amsler et al 20081 ') . Here we consider fermion mass is 
938.3MeV/c^ with spin 1/2, and boson mass is 91.19GeV/c^ 
with spin 1. 

For fermions, Fig.l shows the plot for the parameter 
S in equation (|28p in function of the temperature and the 
chemical potential. The evaluation of the expression for S in 
this case reveals an inverse dependence in square tempera- 
ture and a constant of order 10 to 27, i.e., S ^ W T~ , 
which is the dominant term. In this case, we can see that S 
presents large values for the temperature range. The temper- 
ature range has been chosen to vary from T ~ to T ~ 3K, 
which is the range where significant variations have been 
observed. The chemical potential has been chosen to vary 
between 10~^^ and 10~^^ since in this range the integrals 
present some variation. The statistical mechanics integrals, 
eqs.(13}-(|31), are theoretically defined from zero to infinity, 
but for numerical purposes we used Simpson integration 
with X variating from 0.1 to infinity to avoid a singularity 
in MAPLE during boson computations. 

For bosons. Fig. 2 show that S there is a similar be- 
haviour. Meanwhile, notice that for low temperatures the 
function increases to values that are of order 10^^; this is 
because the expression for S in this case is of order 10^^ T"^. 
One have to point out that despite fermions and bosons 
cases present similar displays, the values on 5 axis are very 
different in order: boson gas can reach higher values, since 
Z° mass is larger than the neutron one. 

If one uses the so expected Higgs boson {H'^, presently 
under search in CERN's LHC), wit h expected mass above 
114GeV (LEP) JAmsler et al 20081 ') . the discrepancies are 
greater. 

The behaviour for the parameter po for fermions is 
shown in fig. 3, and for bosons in fig.4. Aside the integrals in- 
volved, the parameter po scales like T^'^ ^ growing with tem- 
perature. (Observe that the two displays present the axis for 
T and ^ in different direction if compared to fig.l and fig. 2). 

In the range of temperature selected for these sample 




Figure 3. Plot of po function for fermions. Notice that in this 
case po < 2.5. 




Figure 4. Plot of po function for bosons. Notice that the values 
are in order 10^. 
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graphics we have obtained nearly vanishing K functions (for 
both fermions and bosons); meanwhile, for high temper- 
atures we have observed a variating K parameter. In the 
range of temperatures where both S and po have significant 
values the parameter K is nearly vanishing, and vice versa. 
This result for K requires further investigation, since K ^ 
is a necessary condition to avoid a vanishing pressure gas 
(this is a work in progress). 

Anyhow, for any of the numerical ranges tested we ob- 
served that 5 — 1/n (and po) presented a strong correlation 
with temperature and chemical potential, within the limits 
of the expansion used for the equation of state. 



6 LIMITS AND VALIDITY 

The results obtained in the previous section can rise ques- 
tions concerning the expansion validity for any of the numer- 
ical ranges evaluated, since S and po display huge values. In 
fact, figures [T] to |4] form just part of a whole: in those fig- 
ures only domains where the surface presented pronounced 
variation are presented; temperature and chemical potential 
ranges extend to infinity with no limitations up to now. De- 
spite having achieved the aim of this work, which is to say 
showing there is a connection between the polytropic and 
the isothermal sphere models, we perform a deeper analysis 
of the polynomial expansion (|27|) . 

As explained in section [4l equation (|27[) comes from: 

P = KLl+' + {S+l)p',ip-po) + ^-^^-^pl-\p-pof 



+ '-^^PriP~Pof + - 



which is a Taylor expansion of (|25p around the central den- 
sity p = Po- Two points are important to keep the validity 
for this expansion: terms (p — po)^ — > and 5 — >■ 0. 

The first condition is a necessary one otherwise parts of 
the third term in eg. pop mix up with those from the pre- 
vious term. Moreover, the first condition can be interpreted 
as p{r) ~ Po, which is a reasonable demand when consider- 
ing an isothermal gas sphere during the formation of a star, 
for instance. The second condition is necessary since there 
are terms in 5^ that must be neglected in such eq. (|27p can 
be considered an useful approximation and to guarantee its 
resemblance to ea. (|23|l . 

In figures [T] to m we can see that S displays huge values 
and it is questionable if the comparison that results in the 
present model is valid. But, as remarked above the values of 
T and p, extend to infinity. Thence, we can look for regions 
in the domain T x p. where 5 — >■ 0. 

Since 5 = 2a — 1, the limit (5 — >■ yields a — >■ 1/2, or: 



(«-|) 



0. 



Numerically we defined a region where this limit is 
achieved with a certain precision e: 



< e 



(31) 



and we have used this condition to search only for points 
that are within this precision. If we can find such points 
we perform an existence condition. Indeed we were able to 




Figure 5. Search for points attending ea. ll31l l with e = 0.5. 
Crosses correspond to points where 5 — )■ for the fermions case. 
Logarithm is base 10. 



find such points for both configurations, boson and fermion 
stars. This would be expected since 5 decreases as T and p 
increases. 

In fact, considering fermions to be neutrons, then a in 
ea.(f29]l is given by: 



poo y^^ 1+0A592X 10- ^^ Tx 

TT^t:^ 



OPcrmi = 0.1186 X 10^ 



T'l: 



x^/^^l+OA592xlO-^^Tx 



(32) 



f 0. 
:>[x 



7243xl(>;^fj£ ^ 



(30) which shows that, to obtain the required limit we need 



10^'' 1 
«~ -j;r -^ 2' 

and that yields temperatures such as T --^ lO^^K. Since 
eq.((33} demand a big computational effort we use this ana- 
lytical view to search for points in the vicinity of this temper- 
ature (for fermion configurations the program spent about 
5 hours and for bosons about 18 hours, using standard com- 
puters). In the case of Z° bosons we have obtained: 



QBo 



0.1120 X 10'' 



^oo 

Jo 


y^^l+OA725xlO-^^Tx 


oxp(. 0.7243X10^^^^,^ 1 


r^2 r°° 3:5/2^1 + 0, 4725xlO-l=Ta; 


^ Jo 


oxp(^ 0.7243X10^:^^,,-) 1 



(33) 



and for that the required limit for a is achieved with tem- 
peratures T ~ 10"K. 

Fig. [5l then, shows the points obtained for the fermions 
case. Crosses then represent points attending ea. (|3ip with 
e — 0.5, p variating from 1 x lO"'^^ to 2 x 10~^ with step 
2 X 10^^, and temperature variating from 1 x 10^^ to 9 x 10^* 
with step 9 x 10^^ , giving ten thousand points to be tested. 
Hopefully, we obtained some points with e $; 0.5 that are 
those for which the model can serve to the desired purpose 
since 5 is small. Notice that we use a log-log plot due to the 
great discrepancy for the values. 

Fig. |6] shows the points obtained for the bosons case. 
Circles represent points attending ea. (|3ip with e — 0.5, p 
variating from 1 x 10~^^ to 2 x 10~^ with step 2 x 10"^, but 
temperature variating from 1 x lO^'* to 9 x 10^® with step 
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eq.©, numerical values for p (and po) are in lack of a c^ fac- 
tor. Observe that this gives the correct units for mass den- 
sity: notice that p has the same units as 5*2, and once 52 is 
expressed in MKS-units kgm~^ s~^ the adequate correction 
in units of c^ yields kgm"'' as expected. Central densities 
on the fourth row are computed directly using ea. (|28|l and 
are not in correct dimensions. 

From Table [T] one can see that the sample values so 
obtained (last row) give rise to densities that are in ac- 
cordance to those expected in reference texts. In fermion 
stars densities central densities are of order 10'^°kg /m 



IZeldovich fc Novikov WrUt . iRufHni fc Bonazzola 1969 ). In 
boson stars central densitie s are given by po ^ vn? /A -kG 
and are of order 10^**kg/m^ (|Gleiser 19891 : iGleiser 198a ). 



Figure 6. Search for points attending ea. ll31l l with e = 0.5. Cir- 
cles correspond to points where 5 — > for the bosons case. Loga- 
rithm is base 10. 
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Figure 7. Merged plot of both cases, while searching for points 
attending eq. JSTJ with e = 0.5. Spits at the top correspond to 
bosons and spits at the bottom correspond to fermions. 



9 X 10^*, giving again ten thousand points to be tested. For 
bosons we also observed the existence of points for which 
the model apply within the precision e. 

Fig.[7]shows both cases as seen from above. One can see 
some clouds of points (that due to the softness of the func- 
tion grouped like 'spits'). The number of points obtained for 
fermionic configurations was greater than the bosonic ones, 
in the limits tested. From this plot one can see that there 
is a difference at the characteristic temperatures for bosons 
and fermions, and we can define different 'regimes' where 
the model approaches correctly polytropes to isothermal gas 
spheres. 

Once one can rely on the existence of regions where 5 — >■ 
there remains the following question: are corresponding 
values of central density po comparable to those obtained 
in the literature? A rapid test can be made by taking two 
sample points in the clouds of the figure [7] for each case. 
We chose log(p) = —8, and temperatures log(r) = 14.5 for 
fermions and log(T) = 16.5 for bosons. 

In Table [T] numerical values computed for mass density 
are directly listed in the fourth row, for the described sample 
cases. But, during the treatment of relativistic stars, since 



7 DISCUSSION 

In this paper we show that there is a relation between the 
statistical mechanics of a self-gravitating gas sphere model 
and the polytropic model; for short, we propose the termi- 
nology polytropic isothermal gas spheres. We start with the 
statist ical mechanics treatment given for fermion and boson 
gases ijlngrosso fc Rufiini 1972h and find an equivalent equa- 
tion of state that is compared with the polytropic model. 
We observe that this relation is similar to isothermal gas 
spheres, i.e., a polytrope with 7 = 1. The main result starts 
in section |4l where a polynomial expansion is performed and 
in first order in p one can find the suggested correspondence 
for the polytropic model parameters: K, po and 5 = 1/n. 
To illustrate our results we use neutrons as fermions and 
Z*^ bosons. The graphics for 5 for fermions and bosons are 
similar in shape, but very different in values range, since a 
values show up remarkable differences due to the masses of 
the particles. 

Limits of this approach are the following. The main re- 
sult on section |4] is obtained only up to first order. But, we 
can see this matches the models in use for isothermal gas 
spheres. Thence, it is not expected remarkable difi'erences if 
one try a larger polynomial expansion. 

We consider only non-charged particles to avoid further 
interactions. Probably, new physics can be extracted consid- 
ering charges particles like protons, and W bosons or H 
Higgs bosons. 

One advantage in this picture is that eas. (|28|l are de- 
scribing the parameters dependence in temperature and 
chemical potential, and there is no need to separate 
the cases for low tem perature and high temperature 
(jlngrosso fc Ruffini 19721 ). Initially we focused at low tem- 
perature range, closer to the (environment) cosmic back- 
ground temperature, since the considered compact objects 
like neutron stars and boson stars are theoretically shaped 
using zero temperature. 

In order to preserve model validity we performed an 
analysis searching for regions (temperature versus chemical 
potential) where the expansion is adequate. We found those 
regions and corresponding central densities are similar to 
those known for fermion and boson stars. Future work could 
use the set of equations similar to (|28p , with a higher order 
approximation in the expansion, to detail the behaviour of 
the parameters for larger ranges of temperature and chemi- 
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Table 1. Central densities pq: sample cases for boson and fermion numerically obtained by the present 
model compared to those cited in the literature (section |6ll . Quantities dimensions are: fi in Joules, T in 
Kelvins, and po in kg/m^ (fifth row). Attention: there are differences between those po units. Refer to the 
text for further information. 



Case 



PO PO 

(Correct Dimensions) 



Fermions 1x10"* 3.16 X IQi-* 6.5195x10^2 5.87 xlO^'' 



Bosons 1 X 10"* 3.16 x IQi'' 3.7262 x 10^ 



3.35 xlQ-* 



cal potential which could be useful for a more expanded set 
of stellar object s and structure s. 

Merafina iJMerafina 199d ) points out that semi- 
degenerate configurations have had a quite satisfactory 
treatment inside the bulge of the star, but for the outer 
shell, where there is no border, the thermodynamical quan- 
tities expansions fail. This problem can possibly be circum- 
vented if one uses the polytropic approach to describe the 
fermionic stars (neutrons stars and white dwarf stars) and 
boson stars. 

The overall conclusion is that there is a strong cor- 
respondence between polytropic and gas sphere thermal 
models, with an evident correlation within their param- 
eters, e.g., the polytropic index n, and possible temper- 
atures, composition and chemical potentials. This is ex- 
actl y what is expected f r om the theory of stellar interi- 
ors l|Chandrasekhar 19391 : IZeldovich fc Novikov 19711 ) . but 
has never been reported before. Hence, in this paper 
we present an alternative to link parameters, and for 
studying their significant ranges in function of tempera- 
ture and chemical potential. In special , the analysis of 
n — >■ (XJ (JNataraian fc Linden-Bell 19 97). and the envis- 
age of some possible implications for cloud condensation 
l|Honda fc Honda 2003i) is of potential interest. 
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